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Density scaling has a rich history in density functional theory, providing exact conditions for 
use in the construction of ever more accurate approximations to the unknown exchange-correlation 
functional. We define a conjugate potential scaling for each density scaling. This provides exact 
relations on various important density functionals, in particular, relating functionals evaluated on 
exact densities of different potentials. We illustrate these conditions on the two- and four-electron 
ion series. 



I. INTRODUCTION 

Density functional theory (DFT) has become the 
method of choice for many electronic structure calcula- 
tions in both computational chemistry and condensed- 
matter physics [1]. It balances the competing demands of 
accuracy and computational efficiency. The foundation 
of DFT is the Hohenberg-Kohn theorem[2] which states 
that the total ground-state energy of a system can be 
written as a functional of the ground-state density. This 
allows the construction of the Kohn-Sham[3] system of 
non-interacting fermions in a unique one-body potential, 
which is defined to yield the same density as the real 
interacting system. Since the ground-state energy is a 
density functional, this then gives the ground-state en- 
ergy. However, this functional is not known exactly and 
in practice, the unknown exchange-correlation (XC) con- 
tribution must be approximated. 

The plethora of XC functionals in current use in DFT 
is a symptom of the lack of a systematic method for func- 
tional development [4] . In this regard, the more exact con- 
straints that can be discovered, the more tools available 
for developers to use in their quest for more accurate ap- 
proximations. Any useful approximation in DFT should 
try to satisfy these constraints, or be tested against them. 

Density scaling has been a particularly useful tool for 
the analysis and development of DFT. A singular exam- 
ple is uniform coordinate scaling [5], where the coordi- 
nates of a given density are linearly scaled, but normal- 
ization is preserved. This has led to fundamental ex- 
act conditions on the exchange-correlation (XC) energy 
functional [5-8]. For example, the form of the local ap- 
proximation to the exchange energy can be deduced from 
this scaling. The adiabatic-connection formulation[9-12], 
much studied and used in DFT development, is essen- 
tially an integral over the uniform coordinate scaling 
parameter[5, 13, 14]. Here, the electron-electron interac- 
tion is scaled by a constant while the density is kept fixed, 
linking the non-interacting Kohn-Sham and the fully in- 
teracting systems, and leads to many more conditions. 
For example, the adiabatic connection formula is behind 
rationalizing the hybrid approach[15-18]. 



Recently, a different form of density scaling was used 
in the development of the PBEsol functional[19]. Here, 
both the coordinate and the particle number are scaled, 
leading to new insights into the XC functional. We refer 
to this as charge-neutral scaling [20], as it is equivalent to 
simultaneously changing the charges on atoms and the 
number of electrons, so as to keep overall neutrality. 

In this paper, we extend the use of density scaling as a 
tool in DFT. Most importantly we introduce the concept 
that any form of density scaling defines a related form 
of potential scaling. This leads to more exact conditions 
on the various DFT quantities as functionals of densities 
of different particle number. Yang and others [21] have 
emphasized the duality of the potential with the density, 
but have not related scaling of one to the other. 



II. BACKGROUND 

In DFT, the total ground-state energy for electrons in 
a given external potential, v{r) is given as a functional 
of the ground-state density n{r), 



E^[n] = F[n] + V[r 



(1) 



where the external potential energy as a functional of the 
density is 



V[n\ 



d^r n(r)u(r). 



(2) 



and F[n\ is the universal functional which may be defined 
via a constrainted search[22] over all wavefunctions ^' 
that yield density n(r) 



F[n] 



min < *|T + K:c|* >, 



= T[n] + V,M- 



(3) 
(4) 



In the Kohn-Sham method, a reference system of non- 
interacting electrons with the same density is solved, so 
it is useful to write 



F[7i]=T^[n] + U[n]+E^o[n]. 



(5) 



where Ts is the kinetic energy of non-interacting electrons 
of density n(r), U is the Hartree energy and i^xc is the 
exchange-correlation energy, defined by Eq. (5). The 
Hartree energy (sometimes denoted by E-n) is defined as 



U[n] ^^jdPrjd" 



n{v)n{v') 



(6) 



In Thomas-Fermi (TF) theory [23, 24], the universal fimc- 
tional is approximated as 



where T^"' is the local kinetic energy functional, 
with As - (3/10)(37r2)2/3. 



(7) 



(8) 




FIG. 1: The exact radial densities of Beryllium (solid line) [26], 
and of the CN scaled (with C, =2) Helium (dashed line) [27]. 



III. POTENTIAL SCALING 

Consider a density n(r) that is the ground-state den- 
sity of some interacting problem with potential w(r). 
Now, introduce some positive parameter, < 7 < cxd, 
which produces a family of densities, ^^^(r), with 7 de- 
fined so that 7 — > 00 corresponds to the high-density 
limit. A simple example is the uniform coordinate scal- 
ing of Levy and Perdew[5]: 

n^(r) = 7'^n(7r), < 7 < cxd, (9) 

where the prefactor was chosen to keep the density nor- 
malizd. For example, under uniform coordinate scal- 
ing with 7 > 1, the density of He is squeezed into a 
smaller volume, and looks like a distorted version of a 
two-electron ion[25]. This scaling has become a mainstay 
of DFT and leads to many important results. Most im- 
portantly, when particles interact, the coordinate-scaled 
wavefunction is not the ground-state wavefunction of the 
scaled density. Considering such a wavefunction as a trial 
state in the Rayleigh-Ritz principle yields useful inequal- 
ities for the various density functionals[5]: 

T[n^]<l^T[n], 7 > 1, (10) 

Kc[%] >7Kc[n], 7>1, (11) 

and a similar condition applies for the correlation energy 
Ec[n] itself. 

A second example that we focus on here is what we 
call charge-neutral (CN) scaling, in which 

riQir) = C^n(C^/^r), < C < 00 (12) 

and so Nc^ — QN . We use C as the scaling parameter 
to distinguish from coordinate scaling. This choice both 
scales the coordinate and changes the particle number. 



For Coulomb-interacting matter, this ensures neutrality 
as a function of ^. For example, for single atoms, it sim- 
ply implies Zq — QZ and the atom remains neutral. Lieb 
and Simon[28] showed that Thomas-Fermi (TF) theory 
becomes exact for neutral atoms as C ^ 00, and Lieb [29] 
later generalized the proof to all Coulomb-interacting 
matter. In Fig 1, we illustrate this scaling on the He 
atom density. 

In both coordinate and CN scaling, as the scaling pa- 
rameter is taken to 00, the solution simplifies. Under 
uniform coordinate scaling to the high-density limit, the 
system becomes effectively non-interacting. Under CN 
scaling to the high-density limit, Thomas- Fermi theory 
becomes relatively exact. In either case, we can ask how 
the potential changes when the density is scaled. We de- 
fine this as the potential scaling conjugate to the given 
density scaling, but consider it for all values of the scaling 
parameter, not just in the high-density limit. 

Under coordinate scaling, in the large 7 limit, 

v''{r)=-i^v{jv). (13) 

We therefore define our potential scaling by this equation, 
applied for all 7. We use a superscript to indicate that 
the potential has been scaled, not the density. This is 
simply how the external potential would change when the 
density is scaled, if the particles were non-interacting par- 
ticles. For example, for a neutral atom, this changes the 
nuclear charge by 7, keeping the particle number fixed. 
As 7 — > 00, the repulsion between electrons becomes neg- 
ligible relative to the nuclear attraction, and the density 
becomes that of the non- interacting limit, scaled by 7. 

Similarly, under CN scaling with <^ — > 00, the TF equa- 
tions become relatively exact [30], and 

v^{y) = C'*/^w(C^/3r), Nq = CN. (14) 

Again, the conjugate potential scaling is defined by this, 
applied to all values of Analogously, if self-consistent 
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TF theory were exact, this is how the potential would 
scale for any ( as the density is scaled. 

Although chosen to match the corresponding density 
scaling in the high-density or high-potential limit, these 
potential scalings can be applied for any values of their 
scaling parameter. Since scaling the potential is much 
more common in quantum problems than scaling the den- 
sity, often solutions are known or can be accurately cal- 
culated for different scalings of the potential, but not 
of the density. In this paper, we find relations and in- 
equalities between such solutions that complement the 
ground-breaking results of the previous generation [5]. 



IV. UNIFORM COORDINATE SCALING 

In the old work [5], Levy and Perdew compared two dif- 
ferent wavefunctions with the same density, whereas we 
compare two different wavefunctions in the same poten- 
tial. To do this, begin from a given potential v{r) with 
ground-state density n(r). Define n''{r) as the ground- 
state density of ^^'''(r), given by Eq. (13). Then rt^''y^(r) is 
a useful trial density for the original problem. It is found 
by first scaling the potential, solving the problem, and 
then scaling backwards to the original problem. (In Fig. 
1, the dashed line corresponds n^(r) for the He density, 
with ( = 2.) This is exactly what was done (but with an 
approximate scale factor) in Ref. 25. 

If n^^^(r) is used as trial density for v{t), the varia- 
tional principle states that 

F[nJ/^] + 7-V^[n^] > F[n] + j-^V^in^], (15) 

which may be rearranged as 



FK, J - F[n] > j-'iV-'in,] - F^[n^]). 



(16) 



Conversely, n^(r) may be used as a trial density for v^{r), 
yielding 



(17) 



which can also be rearranged as 

F[nT] - F[n^] < V^[n^] - V'^[n'']. (18) 

Combining the two inequalities yields a constraint on the 
universal functional F[n]: 



F[n^] 



(19) 



7 7 
which may be written in a concise form, with A = 1 /•y, 

AF^[nl^^]>AF^[n], (20) 

where 

AF^[n]=F[n]-X'F[n,/x]. (21) 



Now, F[n\ is typically dominated by the kinetic energy 
contribution, but this can be removed, because Tg [n^] = 
^'^Ts[n]. Thus 



(22) 



where Fhxc = U + Exc- This tells us that if we begin 
from, e.g., the lowest value of Z that binds a given N 
electrons, then AE^y.^[n\^^] is an increasing function of 
A. 

Simple results can be extracted from this very general 

formula by taking 7 to be very large. This makes n'^(r) an 
essentially non-interacting density, because the external 
potential dominates. Thus 



n]/^{r) nNi(r), A 



0, 



(23) 



where nNi(r) is the density of the system with only an 
infinitesimal electron-electron repulsion. But AF^xcN 
also simplifies as A — > 0, because all terms scale less than 
quadratically. Thus 

AF^xc M ^ i^Hxc [n], A ^ 0, (24) 
yielding the universal result that 

Fexcf^Ni] > -EhxcN; (25) 

applying to all potentials. For 7 < 1, Eq. (20) is less use- 
ful, as most systems of interest lose an electron when the 
external potential becomes too small. To further simplify 
Eq. (22), we note that both the Hartree and exchange 
energies scale linearly with 7, i.e., 



Eaxirij] = 7£'Hx[n], 

so that 

A£;^,[n] = (l-A)F„x[n]. 
Inserted into Eq. (22), we find 



(26) 



(27) 



F„x [nl] + A'E^ [nl] > E^, [n] + A'E^ [n] , (28) 



where 



A'E^[n] = AE^[n]/{l-X). 



(29) 



The simplest way to test this result is by doing a Kohn- 

Sham calculation without any correlation (such as oep 
exact exchange). Then the correlation contributions van- 
ish on both sides of Eq. (28), and so 



E„jc[n] < £'hxK]/7 < ^^Hx[nNi] 



(30) 



This simplifies even further for the special case of two 
electrons in a spin singlet, where Ex[n] = -~U[n]/2, so 
the inequality becomes a bound on the Hartree energy 
alone: 



U[n] < U[n'^]/j < U[n^ 



(31) 
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TABLE I: The Hartree energies, U, for the hehum iso- 
electronic series as calculated with the oep exact-exchange 
method as implemented in the OPMKS code[32]. We also 
demonstrate how, for two values of atomic number Z' , the in- 
equalities of Eq. (31) with 7 = Z' /Z, are satisfied. Note that 
if 7 < 1, the inequality is reversed. The values for bordering 
values of Z bracket the value of U at atomic number Z and 
these bounds become tighter as Z' increases. 



Z 



1 


0.790970 


3.163880 


15.819400 


2 


2.051538 


4.103076 


20.515380 


3 


3.303373 


4.404497 


22.022487 


4 


4.554137 


4.554137 


22.770685 


6 


7.054819 


4.703213 


23.516063 


10 


12.055315 


4.822126 


24.110630 


20 


24.555661 


4.911132 


24.555661 



TABLE IL Hartree-exchange energies for the beryllium iso- 
electronic series. Values were also calculated with the 
OPMKS code with oep exact-exchange. Also shown are 
two examples of the inequalities of Eq. (30), again using 
7 = Z' jZ. Although the quantities are more complicate that 
those in Table I, the overall trend is the same. 



Ion 


Z 


Ehx 


Z'=10 


Z'=16 


Be 


4 


4.489776 


11.224440 


17.959104 


B+ 


5 


6.119120 


12.238240 


19.581184 


0"+ 


8 


10.893545 


13.616931 


21.787090 


Ne«+ 


10 


14.051482 


14.051482 


22.482371 


S12+ 


16 


23.498356 


14.686473 


23.498356 


Ca«+ 


20 


29.788628 


14.894314 


23.830902 



In Table I, we analyze the above inequality, Eq. (31), 
while in Fig 2, we plot \J\n?\l^ as a function of 7 for 
exact-exchange calculations of the two-electron ion se- 
ries, beginning with H~. Indeed, the function increases 
toward the Bohr atom limit of 5/4, found by inserting 
a doubly-occupied Is Hydrogen atom orbital into the 
Hartree energy. 

To test the exchange contribution in a non-trivial way, 
i.e., Eq. (28), we repeated the calculations for the 
four-electron ion series, this time beginning from Be. 
Again the inequality is satisfied, and the limiting value 
is found by evaluating the Hartree and exchange ener- 
gies of doubly-occupied Is and 2s Hydrogenic orbitals, as 
calculated in Appendix B. These values are reported in 
Table II and plotted in Fig 3. 

Lastly, we can even include extremely accurate es- 
timates of the correlation contributions for the two- 
electron series. We work from the data in Table I of Ref. 
[31]. Since the two-electron ions are generally weakly cor- 
related, one can approximate the scaling of their correla- 
tion energies with a Taylor-series around the high-density 
limit: 

Eo{n\^E'^o\n\+\E^o\n\ (32) 
where E^if' \n\ are scale-invariant functionals. Since = 




5 10 15 

Z' 



FIG. 2: Using the Hartree energies from Table I, Eq. (31) is 
illustrated for 7 = Z' jZ and Z = \. The trend is identical to 
that seen in Table I, however it is clear that the value is ap- 
proaching it's asymptote, 5/4. This is the Hartree energy for 
density consisting of the doubly occupied hydrogen Is orbital. 




10 20 

Z' 



FIG. 3: The Hartree-exchange energies reported in Table 
II are used to illustrate the inequalities of Eq. (30) with 
7 = Z' jZ and Z = 4. Compared to Fig. 2, the value of 
£'hx[m^]/7 is not as fully converged to its asymptote, however 
the maximum value of 7 is 4 times smaller. The asymptotic 
value for this case is 586373/93312 = 6.284, which is found 
by doubly occupying both Is and 2s hydrogenic orbitals with 
Z = A: and calculating Hartree and exchange energies. 



—Eq, + 9i?c ["71/97(7 — 1)[5], and Tc is reported in their 
table, one can solve for these two coefficients. This yields 
a value of -47.6 mH for E^q^ for He, in excellent agreement 
with the value of 47.9 estimated in Ref. [25], and predicts 
a value of -56.1 mH for H~. Using this approximate 
scaling, we can insert all terms into Eq. (22) explicitly 



U Z'=4 Z'=20 1 
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and find their behavior. The numerical corrections to our 
previous results are negligible. 



V. CHARGE-NEUTRAL SCALING 

In this section, we repeat all the logic of the previous 
section, but apply it now to CN scaling. After repeating 
similar steps (given in Appendix A), we arrive at the 
general result: 



where 



AF"[n] = F[n] - a'^^^F[ni/„] 



(33) 



(34) 



and a = Just as we did for coordinate scaling, we 
can refine our inequality substantially. By construction, 
AF"[n] = for F"'"^[n], so we define the useful func- 
tional: 



F^'[n]=F[n]-F''[n 



(35) 



as the #07>Thomas- Fermi contribution to F[n]. Our in- 
equality then reads: 



AF^^«[ny«] > AF^^«[n] 



(36) 



where 



AF^^«[n] = F^^[n] - a^/=^F^^[ni/„] . (37) 

We find an interesting result in the limit a 0, if we 

make the reasonable assumption that all non- Thomas- 
Fermi contributions scale less strongly than C,'^^^ : 



(38) 



as TF becomes relatively exact in the high This in- 
equality is fiendishly hard to test, even in the large ( 
limit. Consider, e.g., the He atom. The corresponding 
TF density is well-known[33] but we would have to eval- 
uate the exact interacting functional on it to find the 
non-TF contribution. All the above results also apply 
directly to non-interacting electrons in a potential, such 
as the Bohr atom[34], with F replaced by Tg, and the TF 
contributions calculated with no Hartree term. But the 
same difficulties remain. 

There is one case where we know enough already to 
test. For the hydrogen atom (or any one-electron sys- 
tem), F = T only, and is given by the von Weizacker 
functional. The TF density (with or without interaction) 
is well-known and singular at the origin, making the von 
Weizacker energy diverge. Thus, the formula is satisfied, 
but not very informative. 

Lastly, we consider Thomas-Fermi-Dirac-Weizsaker 
theory [35] (TFDW). Here we add to TF the local ex- 
change 



I' 



where = — (3/4)(3/7r)^/^, and the next order gradient 
correction to the kinetic energy, 

Both these terms scale the same way under CN density 
scaling, i.e.. 



(r), 



(39) 



F(2)[nd = T(2)K] K] = C^/=^(T(2)[n] +^^^N). 

(41) 

Then can write the inequality as 

F(2)[„C] > ^5/3^(2) j^j^ (42) 

where n(r) has been evaluated self-consistently within 
TFDW and C > 1. Thus 

[n^F] > [n«]/C5/3 > F(2) [n], (43) 

where n^^ (r) is the Thomas- Fermi solution for the same 
potential as for n(r). 



VI. CONCLUSION 

Potential scaling, conjugate to a given density scaling, 
promises to be a useful tool in density functional the- 
ory. It leads to many exact conditions that can be used 
in functional construction. We have applied it to two 
distinct types of scaling: uniform coordinate scaling and 
charge neutral scaling. In both cases, we have found sev- 
eral interesting bounds. Uniform coordinate scaling was 
useful for analyzing Kohn-Sham DFT, leading to inequal- 
ities involving the only unknown in DFT, the exchange- 
correlation functional. The limit of this inequality in- 
volves evaluating the Hartree-exchange-correlation en- 
ergy of the density of non-interacting fermions in the 
external potential. This connection between the interact- 
ing and non-interacting systems resonates with standard 
approaches in many-body perturbation theory. We il- 
lustrate the bounds on the Hartree-exchange energy this 
inequality provides by performing OEP exact exchange 
calculations on helium and beryllium, showing the ap- 
proach to their asymptote. On the other hand, charge- 
neutral scaling provides inequalities involving Thomas- 
Fermi quantities. The Thomas-Fermi approximation be- 
comes relatively exact for all electronic systems[28, 29] 
and these relations link the corrections to Thomas- Fermi 
with the true system, including those of TFDW theory. 
However evaluating the TF density within these theories 
often leads to divergences [36, 37]. It is not clear how 
important these relations will be for functional develop- 
ment, but those derived from uniform coordinate scal- 
ing are siblings to those derived by Levy and Perdew[5], 
which proved very useful in constraining approximations 
in DFT. 

We thank Eberhard Engel for use of the OPMKS code 
and also thank Cyrus Umrigar for providing exact den- 
sities. Wc acknowledge support from NSF under grant 
CHE-0809859. 
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APPENDIX A: CHARGE-NEUTRAL SCALING 
INEQUALITY 



The Hartree energy is defined by Eq. (6), however in the 
special case of spherical densities it may be written as: 



We follow the steps in deriving Eq. (20) but applied 
to the charge neutral scaling defined in Eq. (14). Taking 
nj^^(r) as a trial density for the v{r) system, then the 
variational principle states: 



F[n'l^^] + V[ni/^]>F[n]+V[n] 



(Al) 



where n^(r) is the exact density for the scaled potential 
v^{r). Conversely, use nf(r) as a trial density for the 
v^{r) system: 



(A2) 



where we have used = ^''/"^^[n]. Combining these 

inequalities gives: 



i^K] - F[n'^] 



> F[n] - F[n{.A 



(7/3 

which may be written as 
with 

AF^[n]=F[n]-^. 



(A3) 



(A4) 



(A5) 



APPENDIX B: EXCHANGE ENERGY FOR 
NON-INTERACTING BERYLLIUM 

The limit of the inequality, Eq. (28), is the Hartree- 

exchangc functional evaluated on the density of the corre- 
sponding non-interacting system. Since the g.s. orbitals 
which sum to this density are known analytically (they 
are simply hydrogenic orbitals) , we may calculate the ex- 
act Hartree-exchange value. 

Written in spherical coordinates, r = |r|, the Is and 
2s hydrogenic orbitals are: 



Mr) = — 



-Zr 



(Bl) 



1/2 



= W=l32;^j (2-Zr)e--/^ (B2) 

For beryllium, both these orbitals are doubly occupied, 
giving the total density as 



U[n] = \ I dr {f[n]{r)f , (B4) 

where 



POO 

f[n]{r) =An dr' r'n(r') , (B5) 

J r 



and we use square brackets to indicate that it is a func- 
tional of the density. The exchange energy for a spin- 
unpolarized system is: 



, </)*(r)</.*(r')^,(r)<^,(r') 



r - r' 



(B6) 

where the factor 2 is due to spin, and the sum is over 
occupied orbitals only, in this case Is and 2s. If we define 
a new quantity, n(r): 



n{r) = ^i3(r)^2s(r) , 



then we may write 



= -2 {U[n,^] + 2U[n] + f/K]) . 



(B7) 



(B8) 



We can use Eq. (B5) and Eq. (B4) for each term sep- 
arately and then combine to find the total exchange en- 
ergy. The answer will be equivalent to solving Eq. (B6) 
using the orbitals, however this method avoids perform- 
ing integrals involving l/|r — r'| and, in this case, are 
easy to solve using integration by parts. The values for 
Hartree, exchange and their sum are: 



, 5 X 23 X 431 49565 „ 
= 7^ = TT^ = 9-561 , 



2634 



59 X 71 X 73 



2736 



5184 



305797 
' 93312 



-3.277 



383 x 1531 586373 
E^^ = = = 6.284 . 



2736 



93312 



n(r) = 2|</.,.(r)|2+2|</>,.(r)|^ 



(B3) 
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